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We consider the totally asymmetric simple exclusion process (TASEP) on a finite lattice with 
open boundaries. We show, using the recursive structure of the Markov matrix that encodes the 
dynamics, that there exist two transfer matrices Tl-i,l and Tl-\,l that intertwine the Markov 
matrices of consecutive system sizes: Tl-i,lMl-i = MlTl-i,l- This semi-conjugation property of 
the dynamics provides an algebraic counterpart for the matrix-product representation of the steady 
state of the process. 



PACS numbers: 05.40.-a, 05.70.Ln, 02.50.-r 

I. INTRODUCTION 

Non-equilibrium statistical mechanics has progressed a lot thanks to the study of interacting particle processes in one 
dimension [l|-|3|- I n this field, the asymmetric simple exclusion process (ASEP) has been playing a paradigmatic role 
with an impressive body of knowledge accumulated during the last twenty years^-lT] . This model can be investigated 
from very different points of view: simple exclusion indeed has a complex story [8| . Two noteworthy techniques are the 
Bethe Ansatz [9l. llOj. useful for spectral properties of the dynamics, and the matrix Ansatz, initiated in in which 
the stationary measure of the model is represented as a matrix-product state. The matrix Ansatz first appeared in the 
study of the totally asymmetric simple exclusion process (TASEP) on a finite discrete lattice with open boundaries, 
as a trick to represent the stationary weights of the configurations. It was observed, empirically, [IH, [l3[ that the 
weights of configurations of sizes L and L — 1 are related through recursion relations that can be suitably encoded in a 
matrix-product state. Subsequently, this idea bloomed into a fruitful and powerful technique that can be summarized 
as follows: given a stochastic model, look for a suitable algebra to represent its steady state. A recent and exhaustive 
review of the matrix-product representation for stochastic nonequilibrium systems can be found in [hlj . 

However, the very fact that the weights of configurations of different sizes can be related through some combinatorial 
identities is deeply puzzling; indeed, two models of different sizes have a priori no relation at all with one another: 
their phase spaces are totally disconnected (the dynamics conserves the size of the system). Moreover, the matrix 
Ansatz does not seem to be logically related to the structure of the Markov operator: the algebra used to represent 
the stationary weights has often to be determined by inspection of simple cases or by analogy with known examples 
[l5l - l20j (the review [3] that describes in detail most of the exact solutions in the field is particularly helpful in this 
respect). Once the algebra is found, the steady state is written as a trace over this algebra and is shown to vanish 
under the action of the Markov matrix. This last step of the proof involves a cancellation mechanism usually requiring 
an auxiliary ('hat') algebra [2~lL l23 |. It has been shown rigorously that most of the models admit formally a matrix 
Ansatz but, unfortunately, the proof is not constructive [23j. 

The aim of the present work is to show for the TASEP with open boundaries that the Markov matrices of two 
consecutive system sizes are semi-conjugate of each other through two transfer matrices. This conjugation property is 
a characteristic of the dynamics and it only relies on the recursive structure of the Markov matrix. A similar property 
has been proved for the multi-species exclusion process on a ring for which the matrix-product representation involves 
complicated tensor products of quadratic algebras [24 - I26I ]. Recently, a conjugation property has been used to derive 
exact results for an annihilation model for which a matrix Ansatz could not be found [27]. We believe that the 
existence of a dynamical conjugation that relates a given model to a simpler one (simpler because it involves a smaller 
number of sites, or of particles, or of types of particles...) is a fundamental property that underlies the solvability of 
many nonequilibrium processes. 

II. RECURSIVE STRUCTURE OF THE TASEP MARKOV MATRIX 

We first recall the dynamical rules of the TASEP in a one-dimensional discrete lattice of size L with open boundaries. 
In this model, particles are injected at site 1 at rate a and rejected from site L at rate /?. Every site can be occupied 
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by at most one particle. Particles in the bulk can hop stochastically with rate 1 from a site to the adjacent site on its 
right if it is vacant. The phase space £Il of the system consists of 2 L different configurations. 

The Markov matrix for a general stochastic system on a Id lattice of size L with open boundaries and nearest- 
neighbor interactions is given by [23j 
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Here, we restrict ourselves to a two-dimensional state space. Thus ft^ and ftA r ) are 2x2 matrices reflecting the 
boundary interactions and ft is the local 4x4 matrix for the bulk. More precisely, we have: 
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Throughout this paper we use 1^ for the 2 d dimensional identity but for d 
that the following recursion for the Markov matrix is satisfied 



M L = 1 <8> M L _i - 1 (8) ft 
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1 the index is supressed. We observe 
l + ft L _i(l,2)cg)l. (5) 



We now explain this formula: a system of size L can be obtained by adding a site with index to the system of size 
L — 1 between the left reservoir and site 1. Then, the Markov matrix Ml for the larger system can be expressed 
in terms of Ml-x- Naively, one would write Ml = 1® Ml-x'- this is correct in the bulk but not in the vicinity of 
the left boundary. Therefore this incorrect formula has to be repaired as follows: (i) particles now enter at site so 
we must subtract the matrix elements that correspond to injecting particles at site 1 which no longer occurs (second 
term on the rhs). (ii) We add a matrix which makes the particles enter at site (third term), (iii) Finally, particle 
hopping from site to site 1 is implemented by the fourth term (which is equal to ft ® Hl-2)- 

For the TASEP, substituting the explicit expressions (jl]) in equation {5), we are led to the following recursive 
structure of the Markov matrix: 
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III. TRANSFER MATRICES FOR THE TASEP 

We shall prove that the Markov matrices Ml and Ml+x for two consecutive system sizes are related by the relation 

f L ,L + xM L = Ml+xTl.l+x- (8) 
This semi-conjugation relation can be illustrated by the following commutative diagram: 

M L 

n L <- n L 



Tl,l+x 



L+X 



^L+X 



Ql+x 



M l+ x 
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We recall that Ql and FIl+x are the phase spaces for the systems of size L and L + 1, respectively. The fact that the 
diagram is commutative means that there are two equivalent paths to evolve from a state Cl in Ql- either one first 
applies the dynamics Ml and then identifies the result to a state in £Il+x via Tl.l+i, ° r one imbeds Cl in Ql+i via 
Tl,l+x and then applies the dynamics Ml+i in the larger phase space. It is important to note that the two transfer 
matrices T and T are different. If they were equal then a whole fraction of the spectrum of Ml would be contained 
in Ml+i [26[ and this is not true for the TASEP with open boundaries as can be verified with simple explicit cases: 
the only common eigenvalue shared by Ml and Ml+i is 0, corresponding to the steady state. The fact that the two 
transfer matrices are different is in constrast with the annihilation model studied in (27j where T = f. 

Nevertheless, the very existence of this commutative diagram expresses the fact that the dynamics for sizes L and 
L + 1 are semi-conjugate to each other. This property also allows us to construct recursively the steady state of a 
system of size L + 1 knowing that of the system of size L. If Ml\vl) = 0, the vector \vl+i) defined as 

\vl+i) = T L x+i\vl) , (9) 

satisfies, using flSJ), 

Ml+x\vl+x) = . (10) 

Hence, if \vl+i) is not the null- vector, it is (but for a normalisation factor) the stationary state of Ml+i- This justifies 
the name 'transfer matrix' for Tl.l+x- We also remark that Tl,l+x has played no role in this construction. One could 
use T as a transfer matrix for the left ground state, but for the Markov matrix it is known that the left ground state 
is the line- vector with all components equal to 1. 

We finally emphasize that the commutative diagram that encodes the semi-conjugation property is an intrinsic 
characteristic of the dynamics and the knowledge of the steady state of the system is not required. The transfer 
matrices can be found for small systems by solving a linear system and their existence relies on the recursive structure 
of the Markov matrix itself. 



IV. PROOF OF THE SEMI-CONJUGATION PROPERTY 

In this section we prove the relation ((HJ by constructing explicitely the transfer matrices T and T. The transfer 
matrix Tl,l+i for L > 1 is given by 

IW = ( «lf +i ) with T^ +1 = ( J" *g l ) and T« := I. (11) 

The matrix Tl,l+x has been constructed to mimic the recursive algorithm provided by the quadratic algebra found 
in (Tlj . We also give an expression for Tl,l+i m terms of an unknown square matrix Rl-i- 

f WI . ( ) with fg +I . ( *_-->_ i(Ml i _ n ) . m 

Using these expressions of T and T, we calculate the left hand side and the right hand side of equation ([8]). Both 
sides are block rectangular matrices of size 4 by 2 with elements given in terms of the matrices at level L — 1. In 
order to satisfy equation ©, the 8 elements of the l.h.s. matrix must be equal to the 8 elements of the r.h.s. matrix. 
Amongst the 8 conditions thus obtained, 6 are tautologically true. The 7th relation is given by 

PL-X-I-L-X,L - ■ ( 13 > 

This equation is easily verified by induction by using the explicit expressions (J7J of pl-x, vl—x and that of T L _ 1 L 
given in equation (jll[) . The 8th relation to be satisfied is given by 

Ti-L-x + T { L 2 } 1L (M L -x + Wi_i) = R L -x [apL-x - (Ml-x - v L -x) [M L -x + u L -x)] ■ (14) 
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This relation can be interpreted as a definition of the unknown matrix Rl-i- In order words, the semi-conjugation 
property (J5|) will be proved, if we are able to construct a family of matrices R that satisfy equation (fl4l for each 
system size. We can rewrite the generic equation for the unknown Rl as 

G L = RlH l , (15) 

where we have defined 

G L = 1l + (M L + oj l ) T L % +1 and H L = ap L - (M L - v L ) (M L + U) L ) = v L M L - M L [M L + u£) , (16) 

the last equality resulting the fact that = apL + vl^l as seen from equation ©. 

Using the finer structure (J6]) of the Markov matrix and the expression for T L L+1 given in equation (jlip . we deduce 
the following recursion for G: 

M L - vl + all M L + + Pl 

' L 



S~i / f±Li "Li I "-"-A, I ^ Li I HL \ fll\ 

G L+ i = ( g ) , (17) 



which can be iterated to obtain a well defined upper triangular matrix for G. Therefore all the Gl's are known. It 
remains to extract Rl from equation (|15|) . It is important to note that there is no need to calculate R explicitly, since 
T plays no role in constructing the stationary state (which is entirely determined by the right transfer matrix T) : the 
only thing to prove is that Rl exists, i.e. that the equation Gl = RlHl has at least one solution. Thus equation (fT5")) 
must satisfy a solvability condition which here amounts to the fact that any vector in the (right) kernel of Hl must 
also belong to the (right) kernel of Gl- i.e., any ket such that Hl\Iil) = must satisfy Gl\}il) — 0. By studying 
explicitely systems of small sizes, we found that the matrix Hl is not invertible and that its kernel consists only of 
a one-dimensional vector-space. We also checked that this one-dimensional vector-space is included in the kernel of 
Gl (which is of dimension 2 L — L — 1). Formal proofs of these facts are given in the appendix. This shows that Rl 
exists and concludes the proof of the semi-conjugation property l|8|). 



V. CONCLUSIONS 



To summarize, we have shown that for the TASEP with open boundaries the Markov matrices M for consecutive 
system sizes are related to each other by a semi-conjugation relation ([8]) via two different transfer matrices T and 
T . An explicit form of T is obtained from the matrix-product ansatz and the existence of the left transfer matrix 
T is proved. This relation between the dynamics corresponding to two consecutive system sizes is a consequence of 
the recursive structure of the Markov matrix. We believe that this correspondance expresses a fundamental property 
which is at the heart of the exact solution of many nonequilibrium models. One advantage of this approach is that 
the existence of a semi-conjugation expressed by equation ([5]) can be investigated on small system sizes by solving 
a set of linear equations, whereas we do not know how to test the existence of a matrix-product representation on 
small system sizes. Above all, we were curious to understand the relation between the recursive structure of the 
Markov matrix and the one implied by the matrix-product representation. It seems to us that the present work gives 
a partial answer to this question. Transfer matrices do appear in other related models such as the multi-species ASEP 
on a ring [25[, the recently studied annihilation model with exclusion (27l | and for discrete-time dynamics such as 
ordered sequential and fully parallel update [30] . It would also be of interest to extend this approach to the partially 
asymmetric exclusion process with open boundaries. 



We thank Arvind Ayyer for useful discussions and S. Mallick for a careful reading of the manuscript. Further MW 
likes to thank Prof. M. Schreckenberg for support. 



Appendix 



The matrix Ansatz for the TASEP: 



A configuration C of the TASEP on a discrete lattice of L sites can be specified by assigning the binary values of L 
local variables Tj with Tj = 1 if site i is occupied and = otherwise. In the long time limit, the TASEP reaches 
a steady state with a non-trivial stationary measure: the stationary probability p(C) of the configuration C can be 
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written as a matrix element over an ordered product of L matrices E and D representing empty and occupied sites 
respectively: 

1 L 

P (C) = —(W\H(r i D+(l-T i )E)\V). (18) 



where Zl is a normalisation factor. As shown in the weights defined in (|T5|) correspond to the steady state 
probabilities if the operators E and D satisfy, along with the left and right boundary vectors (W\, \V), the following 
algebra 

DE = D + E, (19) 
(W\E = a- 1 (W\ (20) 
D\V) =r 1 \V). (21) 

Thanks to these reduction relations, the weight of a configuration of size L can be expressed as a linear combination 
of the weights of some configurations of size L — 1. This matrix-product representation yields exact formulae for the 
currents, the density profiles, the steady correlations, and allows to determine the exact phase diagram of the model 
fill ] . This algebra implies the formula (fTTj) for the transfer matrix: the idea is to always apply rule ([2TH) if possible; 
with second priority apply (|T9l) and if this too is not possible then apply (l2~Tj) . 



Proof that the Kernel of Hl is one dimensional: 

The matrix Hl, of size 2 L , is defined in f| 16[) . It can be considered as the 'unevaluated determinant' of the matrix 
Ml+i, given in equation ©. More precisely, we have: 

det M L +i = det(-H L ) = det H L . (22) 

To write this equation we treat Ml+i as a 2 by 2 block matrix whose elements are themselves matrices of size 2 L and 
we apply the following theorem, proved in [28j: if A, B, C and D are square matrices and if C and D commute, then 

det ( C D ) = det ( AD ~ BG ^ ' ( 23 ^ 
Here C = al^. More generally, using the same theorem, the characteristic polynomial ir(X) of Hl can be written as 



ir(X) = det(Xl L - H L ) = det 



M 



L+l 







(24) 



The kernel of Hl is of dimension 1 if tv(X) is divisible by X but not by X 2 . We show now that for small X, ir(X) 

x f 1 L 



is of order X. We recall that ir(X) is given by the product of the eigenvalues of the matrix 



M, 



L+l 







For X = 0, we know that the spectrum of Ml+i consists of the eigenvalue (with multiplicity 1) and that all the 
other eigenvalues have strictly negative real parts (Perron- Frobenius Theorem, see e.g. (29^)- For X sufficiently small 
the non-zero eigenvalues will remain away from 0; the zero eigenvalue of Ml+i will now be given by Eq(X) that can 
be calculated at first order in X using perturbation theory: 

(0\K\0) , rnlv2 , X . . /0 \ L 



Eo(X) = + 0(X*) = -- (0| I q q j |0) + O(X') , (25) 

where (0| = (1, 1, . . . , 1) and |0) are the left and right eigenvectors of Ml+i associated with the eigenvalue (i..e |0) 
is the steady state probability vector). Both these vectors have strictly positive entries and therefore the dominant 
contribution to Eq(X) is of order X. We have shown that the expansion of the characteristic polynomial tt(X) for 
small X begins with a term of order X with a non- vanishing coefficient: this proves that is a simple eigenvalue of 
H L - 

Explicit construction of the Kernel of Hl ■' 

Let us call {/il(ti, . . . ,tl)} the components of a non-zero vector Hl in the Kernel of Hl- These numbers can be 
obtained in a matrix-product form: 

L 

h L (n, ...,t l ) = (w\Y[ [ n D + (l - n )E] \v) , (26) 

1=1 
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where D, E and \V) are the same as in (fT5H!n]) and (W\ = (W\D. One has to show that Hl acting on Kl gives the 
null-vector. The vector Hl can be written as a column-vector 



/ (W\ES L -x\V) 
\ (W\DS L -x\V) 

where Sl-x stands for all possible strings made of L — 1 symbols D and E. Then, we have 



(27) 



f (W\ES L -x\V) \ f (W\S L -x\V) - a(W\ES L -x\V) \ 
L \ (W\DS L -x\V) ) \ a(W\ES L -x\V) - (W\S L -x\V) J' [ > 

This equation is derived by considering the generic expression Sl-i = D ni E n2 . . . E nk (with n% + . . . + rik = L — 1). 
We then observe, as in [31j . that through the action of Ml all terms in the bulk cancel out (thanks to the algebra 
DE = D + E), all terms on the right boundary cancel out (because of the rule (HU) and only the left boundary 
terms do not simplify because the bra-vector is (W\ instead of (W\. It is important to remark that the precise form 
of (W\ has not played any role in the derivation of equation (f28f and that this relation would remain true if (W\ were 
replaced by an arbitrary bra- vector. 

From this result, we deduce, using equation ([7]), that 



(W\ES L -i\V) 
(W\DS L -i\V) 



0. (29) 



Besides, we have 

(M +„ )( (W\ES L -i\V) \_( {W\S L -x\V) - {W\ES L -x\V) \_( {W\ES L -x\V) \ 

where we have used (W\ — (W\D and DE = D + E in the last equality. Combining equations (|2U| and (|50]) . and 
using equation (Ti7j|) . we conclude that 

f (W\ES L -x\V) \_ , , „ ( (W\ES L -x\V) \ _ ( (W\ES L -x\V) \ _ , , 

The last equality is true because the vector on which Ml acts is precisely the TASEP steady state vector, given by 
the algebra ([TOH2Tj) . 



Proof that the Kernel of Hl is contained in the Kernel of Gl- 
Using the recursive formula (fTTf we can write Gl as follows: 

G L = K L M L + ( q ° Gl i ) with kl = ( J^ 1 I 1 - 1 ) . (32) 

From equation (12811 . we deduce that (klMl)Hl — and this would remain true if {W\ were replaced by any bra-vector. 
The action of Gl on Hl is thus given by 

/ (W\ES L -x\V) \ = [00 \f (W\ES L -x\V) \ ( \ 

V (W\DS L -x\V) J \0 G L -x J \ (W\DS L -x\V) J \ G L -x(W\DS L -x\V) ) ' 1 ' 

This identity is satisfied regardless of the precise form of the bra- vector (W\. Finally, rewriting the string Sl-i = 

n < c L_2 I anc ^ defining (W\ — (W\D, the last equation can be recasted in a form that makes its recursive structure 
Uol-2 ) 

clear 



G L h L 



( ^ 

i ( (W\ES L -2\V) \ ] . (34) 
V L 1 V (W\DS L -2\V) 
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Using iteratively equation we conclude that GlIil = 0, i.e. that the kernel of Hl is included in that of Gl and 
that, therefore, the equation ([15)) is solvable. 
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